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I.  INTRODUCTION 


There  are  many  situations  where  there  is  a  need  to  be  able  to  compute  turbulent  flow 
over  steps:  both  forward  and  backward  facing.  Some  examples  from  the  field  of  projectile 
aerodynamics  include  the  flow  over  rotating  bands,  buttress  threads  and  groves  on  ar¬ 
tillery  shell.  Another  example  is  the  internal  flow  around  injector  plates  and  flame-holders 
in  ramjet  propulsion  systems.  The  use  of  Navier-Stokes  computational  techniques  to  com¬ 
pute  the  flow  around  protuberances,  such  as  rotating  bands  on  an  artillery  shell,  has  been 
attempted  in  which  grid  elements  have  been  “blanked-out”  to  simulate  the  band.1  The 
Baldwin- Lomax2  turbulence  model  was  employed  in  this  application  without  modification. 
Although  the  results  were  encouraging  there  were  a  number  of  questions  concerning  the 
applicability  of  the  turbulence  model  to  local  regions  of  separation  caused  by  protuber¬ 
ances. 

The  following  report  attempts  to  consider  some  of  these  conceptual  difficulties  asso¬ 
ciated  with  the  Baldwin- Lomax  method.  The  first  section  investigates  the  basis  for  the 
method  with  the  aim  of  understanding  its  limitations.  In  the  subsequent  sections,  an  alge¬ 
braic  turbulence  model  is  proposed  for  protuberance  induced  separated  flows  which  is  an 
extension  of  the  Baldwin-Lomax  model.  The  new  model  is  incorporated  into  an  explicit, 
axisvmmetric  Navier-Stokes  solver  developed  by  Patel.3  This  solver  uses  the  MacCor- 
mack  predictor-corrector  method  in  a  zonal  scheme  which  is  particularly  well  suited  to 
computation  of  flows  over  geometries  with  sharp  corners.  The  technique  is  then  applied 
to  a  cone-cylinder-protuberance  configuration  for  which  wind  tunnel  pressure  distribution 
measurements  are  available. 


II.  A  BASIS  FOR  THE  BALDWIN-LOMAX  TURBULENCE 

MODEL 


The  underlying  assumptions  are  that  the  turbulent  eddy  viscosity,  fj,t.  is  a  scalar 
variable  dependent  only  on  local  conditions  (i.e.,  diffusion  and  convection  of  turbulence 
are  not  specifically  taken  into  account).  Given  the  eddy  viscosity,  the  effective  turbulent 
stresses  in  the  mean  flow  equations  can  be  computed  for  local  values  of  the  strain  rate. 
All  these  assumptions  are  common  to  most  isotropic  algebraic  turbulence  models.  In 
addition,  the  Baldwin-Lomax  method  is  a  two  layer  model  composed  of  a  wall  layer  based 
on  Prandtl  s  mixing  length  with  Van  Driest  damping  near  the  wall  and  an  outer  layer 
based  on  a  nearly  constant  eddy  viscosity  in  the  fully  turbulent  ‘wake’  region. 

The  wall  layer,  as  formulated  in  reference  (2),  is  described  by  the  following  formula: 

p\  =  p(2  M-  (1) 
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where: 

P\ 

P 

£ 

kl 

y+ 

A+ 

Ut 

y 


inner  layer  eddy  viscosity  , 
density  , 

Ky[\  -  exp(  — y+/.4+)]  , 

where  the  square  bracket  term  is  Van  Driest's  damping  factor 
magnitude  of  the  vorticitv  , 

Vry/v u.  , 

26.0  (Nondimensional  damping  parameter=  uTA/v), 

shear  velocity.  \fru/ pu 
distance  normal  to  the  wall  . 

0-4  (von  Karman  constant)  . 


This  is  the  same  as  the  earlier  Prandtl-Van  Driest  mixing  length  approach  except  that  the 
derivative  of  the  velocity  with  respect  to  the  normal  to  the  wall  is  replaced  by  the  vorticity. 

In  the  outer  region  of  a  viscous  layer  two  formulations  are  proposed,  one  for  attached 
boundary  layers  and  the  other  for  wake-like  separated  flows.  The  turbulent  viscosity  is 
taken  to  be  the  smaller  of  the  following  two  values: 

(1.) 


P°  =  KCcrpymFiyrr.h ,  (2) 

where: 

F{y)  =  J/M[l  —  exp(  —  y+//l+ )]  , 

F{yrrt)  =  maximum  of  F(y)  where  ym  is  the  normal  distance  to  the 

point  where  F(y)  is  maximum  . 

K  =  0.0168  . 

Ccp  =  1.6  . 

7  =  [1  +  5. b(Ckieb  y/j/m)6]-1  (Klebanoff  intermittency  factor), 

Ckieb  =  0.3  (note:  original  form  for  7  implies4  Ckub  =  Vm/S). 
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(2.) 


//;  -  KftCcyC^\U2dym/F{ym)h.  (3) 

where: 

Ud  —  \J(y'2  +  l'2  +  U,2)mar  ~  \]iu'2  +  l'2  +  U’2)m,n  . 

Cu =  0.25  . 

Note  that  the  second  term  in  the  I'd  equation  is  zero  for  all  wall  flows  if  the  no  slip 
boundary  condition  is  applied  at  the  wall.  This  is  because  the  velocity  components  are 
squared.  Such  a  form  is  more  appropriate  for  wakes  or  jet  flows  where  the  Ud  is  the 
difference  between  two  positive  velocities.  For  separated  wall  flows  which  involve  some  up¬ 
stream  flow,  the  minimum  velocity  is  usually  considered  negative  relative  to  the  maximum 
velocity.  If  the  velocity  scale  is  interpreted  as  the  difference  between  the  maximum  and 
minimum  velocity,  then  I'd  is  larger  than  the  maximum  velocity. 

Since  these  equations  are  intuitive,  phenomenological  hypotheses,  they  can  not  be  de¬ 
rived  with  rigor.  The  original  Baldwin  and  Lomax2  paper  did  not  include  any  derivations. 
However,  some  aspects  can  be  related  to  well  known  empirical  descriptions  of  turbulent 
shear  flows.  In  the  following  sections  the  relationship  to  turbulent  velocity  profile  charac¬ 
teristics  is  developed  and  discussed  with  the  objective  of  defining  the  applicability  of  the 
method  to  separated  flows  over  step  geometries. 

A  significant  part  of  the  following  development  has  to  do  with  the  evaluation  of  the 
parameters  C'cp.  Cuk  (or  more  correctly  the  product  CcpCu, *)  and  Ckieb-  These  parameters 
are  shown  to  be  related  to  the  parameters  of  the  boundary  layer  profiles  considered.  How¬ 
ever.  it  should  be  kept  in  mind  that  Baldwin  and  Lomax  determined  the  values  for  these 
parameters  by  comparing  numerical  results  with  data  from  a  transonic  shock-separated 
flow  experiment.  Thus,  the  parameters  are  probably  not  universal,  but  depend  on  the 
Mach  number  range,  pressure  gradient  and  possibly  other  conditions  which  have  yet  to  be 
identified. 


1.  ATTACHED  BOUNDARY  LAYER 

A  starting  point  is  to  consider  incompressible,  two-dimensional  boundary  layers.  In 
that  case.  C'lauser5  has  established  that  the  outer  layer  eddy  viscosity  can  be  described 
by  the  formula: 

H°  =  KpUe6 (4) 

where 

61=  [\l-v/V')dy, 

J  O 

IC  =  0.016S. 

This  formula  is  usually  used  to  determine  the  maximum  of  the  eddy  viscosity.  Equation 
(4)  is  then  multiplied  by  the  Ivlebanoff  intermittency  factor  to  account  for  the  decrease  in 
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as  the  inviscid  external  flow  is  approached. 

Coles6'7  has  proposed  the  following  single  empirical  equation  for  the  velocity  profile 
which  unifies  the  wall  and  wake  regions  except  for  the  region  immediately  adjacent  to  the 
wall  (the  laminar  sublayer). 


u  =  lu¬ 


ll 


In  (y+)  +  B  +  2—  sin:  (o.5"| 


(5) 


where: 

B  =  wall  layer  constant,  normally  5.0  to  5.5, 


fl  =  wake  layer  parameter  (function  of  pressure  gradient). 


t  =  boundary  'aver  thickness  . 

The  main  difficulty  with  this  equation  is  in  the  proper  evaluation  of  6  which  is  somewhat 
arbitrary  because  of  the  asymptotic  disappearance  of  viscous  effects  at  large  distances  from 
the  wall.  The  Baldwin-Lomax  method  is  an  attempt  to  circumvent  the  necessity  for  the 
direct  evaluation  of  the  boundary  layer  thickness. 


In  a  wall  bounded  viscous  flow,  sufficiently  far  from  the  separation  point  so  that 
over  most  of  the  boundary  layer  the  Van  Driest  damping  factor  is  unity,  the  vorticity  is 
primarily  du/dy.  Tims,  the  function  F{y)  can  be  formed  as  follows: 


F(y) 


(6) 


where  ?/  =  0 ,b~y/t  .  This  quantity  has  the  dimensions  of  velocity  and  primarily  describes 
the  moment  of  wake  vorticity  with  respect  to  the  surface.  The  maximum  in  F[y)  is  defined 
by  a  unique  value  of  ?/  denoted  by  ?/\  The  value  of  >/’  can  be  determined  from  the  following 
implicit  equation  which  is  obtained  by  setting  dF/dy  =  0. 


sm  >/  cosy 


(i- 


!  sill  r/ 


>/*  -  0. 


Thus  ?/'  is  found  to  be: 


y'  =-  1.014375. 


Since  y  —  y,r<  at  that  point,  t  is  seen  to  be  directly  proportional  to  ym. 


d  =  1.54S  IJm- 


(7) 


This  result  determines  the  constant  in  the  internrittency  function  for  attached  flows  as: 


CkU  6  =  ~  =  0.64G0 
c 


(8) 


Coles  profile  law.  equation  (5).  can  be  integrated  to  provide  a  relationship  between  6'k  and 
6  as: 


=  6 


(i  +  m 


(9) 
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Thus 


1+n  77  yrr 


(10) 


Finally  uT  can  be  eliminated  in  favor  of  F(y,n )  by  using  equation  (6)  so  that  Clauser’s 
wake  eddy  viscosity  becomes  identical  with  equation  (2)  if: 


,  _  '  (1  +  II)  7T 

cp  2  [1  4-  411  (sin  rjm  cos  77*)  r/*]  77* 


(11) 


The  pressure  gradient  parameter.  II.  cannot  be  specified  in  a  general  way.  Coles  estimated 
II  to  be  0.55  for  flat  plates  at  high  Reynolds  number.  For  equilibrium  turbulent  boundary 
layers.  White8  has  correlated  II  in  terms  of  the  external  flow  pressure  gradient  from  a 
number  of  experiments  in  the  following  form: 

n  =  0.8(J  4-  0.5)° 75  ,  (12) 


where 


3  = 


/  dpt 
\TUJ  dx  ' 


and  dpt/dr  —  the  pressure  gradient  at  the  °dge  of  the  boundary  layer.  For  accelerating 
flows,  the  pressure  gradient  and  3  are  negative  with  the  limiting  condition  being  11=0. 
On  the  other  hand,  for  decelerating  flows  with  positive  pressure  gradients,  there  is  no 
limit  to  the  magnitude  of  3  or  II  because  ru.,  the  wall  shear  stress,  approaches  zero  as  the 
flow  approaches  separation.  Thus,  the  limiting  conditions  of  II  for  attached  flows  can  be 
identified  and  the  corresponding  values  of  Ccp  calculated  from  equation  (11). 


n  =  0.  Ccp  = 
n  =  0.55.  ccp  = 
n  =  oc.  ccp  = 


-  =  1.549 
2 

1.2 

77 

8  sin  ?/"  cos  )2 


0.S510  [ 


(13) 


This  is  to  be  compared  to  the  value  of  Ccp  of  l.G  recommended  by  Baldwin  and  Lomax. 


An  analysis  similar  to  the  above  has  bepn  reported  by  Granville.9  He  employed  a 
different  form  for  the  Coles  wake  function  which  was  developed  by  Moses.10  In  addition,  a 
correlation  by  Nash11  was  used  to  connect  the  wake  parameter  17  to  the  pressure  gradient. 
The  consequence  of  Granville's  analysis  is  to  show  that  Ccv  varies  with  II  in  basically  the 
same  manner  as  obtained  here  in  equation  (11).  However,  the  form  of  the  wake  function 
Granville  used  results  in  CkUb  (i.e.  ym/d)  also  being  a  function  of  n  which  is  not  the  case, 
when  the  original  Coles'  wake  function  is  used. 


2.  INCIPIENT  SEPARATED  FLOW 


It  is  assumed  that  the  velocity  profile  for  incipient  separated  flow  can  be  described 
by  the  limiting  form  of  equation  (5)  as  n  becomes  large  and  the  wall  shear  stress,  rw  (i.e., 
ur )  goes  to  zero.  The  logarithmic  and  constant  wall  function  terms  go  to  zero  but  the 
coefficient  of  the  sine  term  has  the  following  limiting  value: 

,  2urU 
inn  - =  L/f. 


— 0 


(14) 


Thus,  the  incipient  profile  can  be  written  as: 


u  =  Ue  sin2  ^ 


Try_\ 
46  )  ' 


where  6  has  been  redefined  in  terms  of  6.  the  half  width  of  the  shear  layer.  In  this  case  the 

y*  ,, 


displacement  thickness  equals  the  half  width,  6.  The  function  F(y)  becomes: 

F(y)  =  y^r  =  2C*e  (sin  r?  cos  r?)  r/. 

F(ym)  corresponds  again  to  y  =  =  1.014375  and 

b  =  *W  («) 

The  Clauser  wake  eddy  viscosity  becomes  the  same  as  equation  (2): 

=  KCeppymF(ym)f. 

where  Ccp  =  0.8510  which  is  identical  to  the  II  — ►  oc  limit  for  the  attached  boundary  layer 

case,  as  computed  earlier  in  equation  (13). 


3.  FULLY  SEPARATED  FLOWS 

Next  consider  a  full.'  •  arated  region  in  which  the  shear  layer  is  displaced  a  distance 
h  away  from  the  wall  ;■  ■  wn  in  the  following  sketch.  The  shear  layer  of  width  2b  is 
now  taken  to  extend  from  th~  minimum  to  the  maximum  velocity.  The  location  of  the 
maximum  in  F{y),  relative  ..e  surface,  is  no  longer  proportional  to  6,  the  half  width  of 
the  shear  layer. 

The  shear  layer  velocity  distribution  is  assumed  of  the  same  form  as  the  incipient 
profile  but  translated  a  distance  h  from  the  wall.  Thus,  the  sine  squared  profile  starts  at 
the  minimum  velocity  location  rather  than  the  wall.  The  shear  layer  profile  is  approximated 
by: 

u  =  umtn  +  ud  sin2  ~  h)]  ; 

46 


6 


for  y  >  h.  where  I'd  =  l  mar  —  t  mm  (£  min 


<  0). 

yi 


The  maximum  velocity  gradient,  du/dy.  occurs  at  y  —  h  -f-  6  so  that. 


b  = 


7T 


t'd 


(16) 


This  relation  can  be  put  into  a  form  consistent  with  that  of  Baldwin  and  Lomax  by  taking 
ym  %  h  +  b  and  F(ym)  «  (h  +  b){du/dy)max-  Therefore,  the  length  scale  for  the  outer  flow 


can  be  written  as: 


,  d  J/m 

0  %  - - - . 

4  F{ym) 


(17) 


Thus  the  turbulent  viscosity  becomes  (see  Equation  (3)): 


K  pC cpCuk 


V2dVm 

F(ym) 


where  CcpCuk  =  0.25*  and  if  Ccp—  0.S510.  then  CV/j  =  0.9229.  This  should  be  compared  to 
the  value  recommended  by  Baldwin  and  Lomax  of  0.25.  Although  a  value  of  Cu.i  =  1.0  is 
currently  in  use  in  some  versions  of  the  implicit  Parabolized  Navier-Stokes  code. 


4.  SUMMARY 

In  the  above  sections,  the  outer  layer  eddy  viscosity  formulas  of  Baldwin  and  Lomax 
have  been  related  to  the  velocity  profile  characteristics  through  the  use  of  Coles  law  of 
the  wall  and  wake.  This  was  done  in  order  to  show  how  the  Baldwin-Lomax  formulas  can 
be  derived  from  existing  semi-empirical  information  and  to  clarify  the  approximations  and 
limitations  of  their  equations.  For  separated  flow,  a  significant  approximation  is  made  in 
going  from  equation  (16)  to  equation  (17).  Thus,  in  the  following  application,  the  outer 
layer  characteristic  length  will  be  computed  based  on  equation  (16).  This  is  consistent 
with  the  spirit  of  the  Baldwin-Lomax  method  but  avoids  some  of  the  limitations. 


III.  MODIFICATION  FOR  STEP  FLOWS 


In  applying  the  above  theory  to  the  flow  over  forward  and  backward  facing  steps 
several  new  problems  arise.  Some  of  these  are  that:  1)  the  F(y )  function  can  have  sev¬ 
eral  maxima  and  minima  and  it  is  difficult  to  numerically  program  to  find  the  correct 
maximum,  2)  wall  damping  effects  must  be  provided  from  more  than  one  nearby  wall, 
3)  the  implementation  of  the  intermit tency  factor  should  account  for  the  displacement  of 
the  shear  layer  away  from  the  wall  by  the  recirculation  zone  and,  4)  the  Van  Driest  wall 
damping  factor  should  not  go  to  zero  at  the  separation  point. 

The  Navier-Stokes  code,  in  which  the  new  turbulence  model  is  implimented,  uses  a 
grid  which  divides  the  flow  region  along  the  body  into  several  zones.  In  two  dimensional 
or  axisymmetric  flows,  the  grid  can  be  designed  so  that  the  step  geometry  is  defined  in 
terms  of  constant  grid  lines.  The  x  coordinate  runs  parallel  to  the  main  flow  direction 
and  y  normal  to  it.  The  code  marches  in  time  and  thus  the  calculation  of  the  turbulent 
viscosity  lags  one  time  step  behind  the  solution  of  the  mean  flow  equations. 

On  walls  with  attached  flow  the  conventional  Baldwin-Lomax  method  is  used.  When 
any  reversed  flow  (i.e..  negative  velocity  relative  to  the  freestream)  is  detected  at  any  given 
station,  the  separated  flow  computation  is  invoked.  Reversed  flow  exists  when  the  velocity 
component  parallel  to  the  local  surface,  at  any  point  on  a  normal  from  the  wall,  is  negative 
relative  to  the  main  flow  direction. 


1.  OUTER  EDDY  VISCOSITY 

The  difficulty  with  multiple  maxima  is  eliminated  by  returning  to  the  basis  for  the 
separated  flow  length  scale.  Instead  of  defining  it  in  terms  of  the  ratio  of  ymUd/F{ym)  (see 
Equation  (3)),  one  half  of  the  ratio  of  the  maximum  velocity  difference  to  the  absolute 
value  of  the  maximum  vorticity  is  used.  The  search  for  the  maximum  vorticity  starts  at 
the  minimum  in  the  velocity.  This  assumes  that  the  vorticity  associated  with  the  free  shear 
layer  is  the  largest  and  the  shear  layer  determines  the  outer  length  scale.  Thus  the  outer 
length  scale  is  calculated  from: 

6  =  0.5-^-,  (IS) 

|u'  (mar 

and  the  outer  eddy  viscosity  is  given  by: 

y°t  =  KCcpUdb (19) 

2.  INNER  LAYER  EDDY  VISCOSITY  FOR  FORWARD  AND  REARWARD 
FACING  SURFACES 

Near  any  lateral  wall,  the  Prandlt-Van  Driest  inner  layer  eddy  viscosity  takes  over 
whenever  it  is  less  than  the  outer  layer  viscosity.  The  surfaces  normal  to  the  main  flow 
direction  also  introduce  a  wall  layer  effect  and  the  corresponding  eddy  viscosity  associated 
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with  such  surfaces  is  calculated  as. 


Tt 


1  —  exp 


2 


(20) 


where  xr.  is  the  normal  distance  for  these  surfaces.  The  eddy  viscosity,  then,  is  the  smaller 
of  three  quantities:  n°.  jj'yt  or  y'It  . 


3.  SEPARATION  POINT 


Computation  of  the  inner  layer  eddy  viscosity  and  determination  of  F(ym )  involves 
an  additional  problem  when  passing  from  the  attached  to  separated  flow  because  the  wall 
shear  stress  goes  to  zero  at  that  point.  The  Van  Driest  damping  parameter,  A,  also  goes 
to  zero  and  the  overall  viscosity  reduces  to  an  unrealistically  low  value  associated  with 
laminar  flow.  This  problem  was  observed  by  Yisbal  and  Knight.12  In  the  present  case,  the 
computed  wall  damping  layer  height,  A.  is  restricted  as  follows: 


The  added  term.  G/6,  remains  finite  when  uT  goes  to  zero.  Thus,  the  maximum  value  of 
A  is  limited  to  1  /G  of  the  total  thickness  of  the  shear  layer,  6.  A  tentative  value  of  three 
has  been  selected  for  G.  The  value  of  6  is  taken  as  the  normal  distance  from  the  wall  to 
where  vt  occurs. 


4.  INTERMITTENCY 


A  problem  arises  in  applying  the  intermittency  factor  because  the  shear  layer  can  be 
displaced  some  distance  away  from  the  wall.  If  ->  is  calculated  as  a  function  of  y/6,  the 
intermittency  is  too  high  at  the  outer  edge  of  the  shear  layer.  This  can  be  corrected  by 
making  =  1.0  for  all  values  of  y  less  than  the  ym  —  b,  where  ym  is  now  the  location  of  the 
maximum  vorticity.  At  larger  values  of  y  we  use: 


=  <  1.0  +  5.5 


(y  -  Vm  +  b ) 
2b 


•]  6  ' 


-1 


y  >  (ym-  b). 


(22) 


Thus  the  intermittency  factor  starts  to  reduce  the  eddy  viscosity  at  the  inner  boundary  of 
the  shear  layer  and  falls  to  approximately  the  same  level  at  the  shear  layer  outer  edge  as 
it  would  at  the  outer  edge  of  an  attached  boundary  layer. 


5.  CONVEX  CORNER 

Another  difficulty  occurs  near  the  top  of  a  forward  or  backward  step  because  it  is  not 
clear  how  to  evaluate  the  inner  layer  eddy  viscosity  in  that  region.  As  previously  described, 
the  inner  layer  vorticity  and  Van  Driest  damping  are  computed  along  a  normal  to  the  local 
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wall.  However,  these  properties  in  the  flow  opposite  a  convex  corner  (such  as  above  and 
ahead  of  the  corner  of  a  forward  facing  step)  can  not  be  reached  by  a  normal  to  either  wall. 
In  order  to  prevent  there  being  a  large  discontinuous  change  in  the  eddy  viscosity  in  this 
region,  the  normal  distance  (xn  in  Equation  (20))  is  taken  as  the  radial  distance  from  the 
corner.  The  shear  velocity  is  allowed  to  change  linearly  with  the  angle  from  one  surface  to 
the  adjacent  one. 


IV.  ROTATING  BAND  APPLICATION 


Experiments  have  been  performed  which  measure  the  pressure  caused  by  supersonic 
flow  over  a  rectangular  band  on  a  cone-cylinder  wind  tunnel  model.13  A  comparison  has 
been  made  between  the  predicted  pressure  distribution  and  these  observations.  Other 
characteristics  of  the  flow  such  as  velocity  profiles  and  contour  plots  of  the  mean  flow 
variables  are  used  to  illustrate  the  quality  of  the  results  of  the  analysis. 

Figure  1  shows  the  grid  used  for  the  computations  and  it  also  indicates  the  experi¬ 
mental  geometry.  The  wind  tunnel  model  consists  of  a  13.1°  half  angle  cone  followed  by 
a  seven  caliber  cylinder  (1  caliber=25.4  mm).  The  Mach  number  was  3.0  and  the  wind 
tunnel  Reynolds  number  was  1.3  x  106.  based  on  body  diameter.  In  the  experiments  a 
rectangular  cross  section  band,  of  thickness  H  =  1  mm  and  width  of  1.27  mm.  was  mov¬ 
able  relative  to  fixed  pressure  taps.  This  allowed  the  pressure  disturbance  caused  by  the 
band  to  be  measured  with  considerable  detail  in  the  nearly  constant  pressure  region  on 
the  cylinder. 

The  computed  pressure  distribution  is  compared  with  the  measurements  in  Figure  2. 
The  agreement  is  very  good  ahead  of  the  band,  near  the  forward  facing  step,  but  un¬ 
derpredicts  the  pressure  level  behind  the  band.  The  minimum  pressure  is  in  fairly  good 
agreement  with  the  data  but  after  the  initial  rise,  the  recovery  toward  the  nearly  ambient 
value  is  computed  to  be  slightly  less  than  that  observed  in  the  experiment. 

Figure  3  shows  the  Mach  number  contours  obtained  by  combining  the  results  from 
the  three  separate  grids. 

Pressure  contours  ahead  of  and  behind  the  band  are  shown  in  Figures  4  and  5.  These 
show  the  compression  and  expansion  waves  generated  by  the  flow  moving  over  the  separated 
zones.  As  expected  the  pressure  on  the  forward  face  of  the  band  reaches  a  maximum  value 
of  three  times  ambient  near  the  top  of  the  band  where  the  shear  layer  flow  reattaches. 
The  mean  pressure  on  the  front  face  of  the  band  is  approximately  15%  greater  than  the 
pressure  on  the  cylinder  just  in  front  of  the  band.  The  pressure  on  the  back  face  of  the 
step  is  nearly  constant;  close  to  the  minimum  pressure.  Except  near  the  leading  edge,  the 
pressures  on  top  of  the  band  are  very  close  to  the  ambient  value. 

The  velocity  vector  plot  of  the  flow  in  the  three  zones  is  shown  in  Figure  6.  This 
plot  clearly  shows  tne  recirculation  regions  near  the  band.  Finally  velocity  profiles  at  two 
stations  are  shown  in  Figures  7  and  8.  Figure  7  shows  the  profile  at  a  point  0.438  step 
heights  ahead  of  the  band.  The  velocity  is  seen  to  fall  very  rapidly  near  y  =  0  to  its 
minimum  value  of  -250  ft/s  and  then  increase  more  slowly.  The  second  of  these  figures 
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shows  results  close  behind  {x/H  =  0.10)  the  band.  Note  the  complicated  variation  of 
velocity  close  to  the  cylinder  wall  followed  by  a  very  thin  shear  layer  coming  off  the  top  of 
the  band.  The  viscous  boundary  layer  is  above  the  thin  shear  layer.  Since  the  turbulent 
viscosity  in  the  inner  layer  is  essentially  proportional  to  du/dy ,  it  exhibits  several  local 
maxima  at  the  inflection  points  and  zero’s  at  the  minimums  and  maximums.  Thus  it  is 
not  surprising  that  the  search  for  the  correct  maximum  in  the  Baldwin-Lomax  function, 
F(y),  is  difficult  to  incorporate  into  a  numerical  code. 

V.  CONCLUDING  REMARKS 


The  Baldwin-Lomax  algebraic  eddy  viscosity  model  has  been  reviewed  and  a  basis  for 
the  separated  flow  formulation  has  been  developed  using  an  assumed  form  for  the  separated 
flow  velocity  profile.  This  analysis  of  the  Baldwin-Lomax  model  has  been  undertaken  in 
order  to  clarify  its  applicability  and  its  limitations  with  regard  to  the  use  in  Navier-Stokes 
computations  of  step  induced  separation.  One  conclusion  which  can  be  drawn  is  that  the 
length  scale  is  associated  with  the  ratio  of  velocity  scale  to  the  maximum  velocity  gradient 
(or  vorticity).  Numerical  computations  can  be  based  on  this  formulation  rather  than  using 
the  Baldwin-Lomax  variables,  F(ym )  and  ym. 

Several  other  modifications  of  the  Baldwin-Lomax  method  have  been  suggested  and 
implemented  for  the  flow  over  forward  and  backward  facing  steps.  These  include:  1) 
providing  for  damping  effects  on  the  forward  and  backward  faces  of  the  step,  2)  accounting 
for  the  displacement  of  the  shear  layer  in  implementation  of  the  intermittency  factor,  3) 
modification  of  the  Van  Driest  damping  function  so  that  it  does  not  become  unrealistic  at 
the  beginning  of  separation  or  reattachment  . 

The  modified  turbulence  model  was  incorporated  into  an  explicit  Navier-Stokes  code 
which  is  designed  to  compute  the  flow  in  zones.  This  feature  makes  the  code  particularly 
suitable  for  application  to  the  flow  over  steps.  The  test  case  which  has  been  computed 
using  this  code  is  that  of  a  cone-cylinder  model  with  a  simulated  rotating  band.  The 
rotating  band  provides  a  rectangular  forward  and  rearward  facing  step  to  the  local  flow. 
The  configuration  has  been  tested  in  a  wind  tunnel  and  detailed  pressure  distributions  are 
available  for  comparison  with  the  computed  results. 

The  results  of  the  comparison  show  that  the  pressure  distribution  ahead  of  the  band  is 
well  predicted.  The  extent  of  the  computed  separated  zone  agrees  with  experiment  as  well 
as  the  form  and  magnitude  of  the  pressure  distribution.  Behind  the  band  the  results  are  less 
satisfactory  in  that  the  pressure  recovery  from  the  minimum  is  too  rapid.  The  minimum 
pressure  and  the  extent  of  the  recovery  region  are  correct  but  the  pressure  at  reattachment 
and  downstream  are  slightly  lower  than  experimentally  observed.  This  result  agrees  with 
other  Navier-Stokes  computations  using  the  un-modified  Baldwin-Lomax  method  with  a 
grid  blanking  technique  to  model  the  rotating  band. 

Additional  testing  of  the  turbulence  model  on  other  flow  situations  will  be  required  in 
order  to  tune  the  empirical  constants  and  investigate  the  range  of  the  models  applicability. 
The  results  shown  here  indicate  that  simple  algebraic  models  can  be  used  in  at  least  some 
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complex  flows. 


Figure  1.  Computational  grid  and  model  geometry  (dimensions  are  in  feet). 


Figure  4.  Pressure  contours  in  separated  flow  region  ahead  of  band 

(dimensions  are  in  feet). 
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Figure  5.  Pressure  contours  in  separated  flow  region  behind  band 

(dimensions  are  in  feet). 
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Figure  7.  Velocity  profile  ahead  of  band  (A  X/H  =  0.438). 
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Figure  8.  Velocity  profile  behind  of  band  (A  X/H  =  0 
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List  of  Symbols 

=  Van  Driest  damping  constant 
=  constant  in  law  of  the  wall 
=  half  width  of  shear  layer 
=  constant  (attached  flow) 

=  constant  in  intermittencv  function 
=  constant  (modifies  Ccr  for  separated  flows) 

=  moment  of  vorticity 

=  displacement  of  shear  layer  from  the  wall 
=  Clauser  constant 
=  mixing  length 

=  Reynolds  number  based  on  body  diameter 
=  pressure 

=  y/(u*  +  v2  +  'r-)max  -  yj(v2  +  v*  +  w2)mtn 
=  x-oomponent  of  edge  velocity 
=  local  x-component  of  velocity 

=  shear  velocity  =  ^Au/pu 
=  local  v-component  of  velocity 
=  local  z- component  of  velocity 
=  distance  along  the  wall 
=  distance  normal  to  the  wall 

=  distance  normal  to  the  wall  to  the  maximum  in  F(v) 


Greek  Symbols 

3  =  pressure  gradient  parameter 

p  =  intermittencv  factor 

6  —  boundary  layer  thickness 

6 1  =  incompressible  displacement  thickness 

;/  =  rry/2d 

k  =  von  Karman  constant 

//  =  dynamic  viscosity 

v  =  kinematic  viscosity 

II  =  Coles'  turbulent  boundary  layer  wake  parameter 

p  =  density 

r  =  shear  stress 

u.‘  =  vorticity 


Superscripts 
i  =  inner  layer 

o  =  outer  layer 

+  =  nondimensional  using  wall  variables 
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Subscripts 
e  : 

min  = 

max  = 

n  '  = 


t 


w 


boundary  layer  edge 

distance  from  the  wall  to  minimum  in  a  variable 

distance  from  the  wall  to  maximum  in  a  variable 

normal 

turbulent 

wall 
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